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ABSTRACT 

Cosmological models with different types of Dark Energy are becoming viable alternatives for 
standard models with the cosmological constant. Yet, such models are more difficult to analyze 
and to simulate. We present analytical approximations and discuss ways of making simulations 
for two families of models, which cover a wide range of possibilities and include models with 
both slow and fast changing ratio w = p/p. More specifically, we give analytical expressions for 
the evolution of the matter density parameter f2 m (z) and the virial density contrast A c at any 
redshift z. The latter is used to identify halos and to find their virial masses. We also provide 
an approximation for the linear growth factor of linear fluctuations between redshift z = 40 and 
z = 0. This is needed to set the normalization of the spectrum of fluctuations. Finally, we discuss 
the expected behavior of the halo mass function and its time evolution. 

Subject headings: methods: analytical, numerical - galaxies: clusters - cosmology: theory - dark energy 



1. Introduction 

Observations of high redshift supernovae (Perl- 
mutter et al. 1999; Riess et al. 1998) as well 
as the analysis of fluctuations of the cosmic mi- 
crowave background combined with data on the 
large-scale structure of galactic distribution (e.g. 
Balbi et al. 2000; Tegmark, Zaldarriaga, & Hamil- 
ton 2001; Netterfield et al. 2002; Pogosian, Bond, 
& Contaldi 2003; Spergel et al 2003) indicate 
that there is a significant component of smooth 
energy with large negative pressure, characterized 
by a parameter w = p/p SS —0.5. This component 
is dubbed dark energy (DE). The nature of DE 
is open for debate with candidates ranging from a 
cosmological constant A to a slowly evolving scalar 
field (f> to even more exotic physics of extra dimen- 
sions (e.g., Dvali, & Turner 2003) 

One of the most appealing ideas for DE is a self- 
interacting scalar field, which evolves with time 
(Ratra & Peebles 1988; Wetterich 1988). We 



call this dynamical Dark Energy. The advantage 
of the dynamical DE models as compared with 
the ACDM models is that DE naturally yield an 
accelerated expansion easing the problem of fine 
tuning. The observational signatures of dynami- 
cal DE should be carefully investigated in order 
to determine which measures can be used to dis- 
criminate ACDM from dynamical DE and among 
different dynamical DE models. In this paper we 
focus on the two most popular variants of dynam- 
ical DE. Ratra & Peebles (1988, RP hereafter) 
studied DE models, which cause a rather slow evo- 
lution of w. Models based on simple potentials 
in supergravity (SUGRA) result in faster evolv- 
ing w (Brax & Martin 1999, 2000). Together RP 
and SUGRA potentials cover a large spectrum of 
evolving w. The potentials are written as 



V{<f>) 



A 4+a 



(RP), 



(1) 



1 



A 4+a 



exp^TrG^ 2 ) (SUGRA). (2) 



Here A is an energy scale, currently set in the 
range 10 2 -10 10 GeV, relevant for fundamental in- 
teraction physics. The potentials depend also on 
the exponent a. Once the parameters A and a 
are assigned, the DE density parameter Qde fol- 
lows. Here, however we prefer to use A and CIde 
as independent parameters. 

Dynamical DE has a kinetic and a potential 
components, reading <ft 2 /2 and V{<j>), respectively. 
Those factors define the energy density pde and 
the pressure pde- In general, the ratio of the pres- 
sure and the density 



Pde = <T/2 - V(4>) 

PDE ft/2 + V{(j)) 



(3) 



changes with time and is typically negative when 
the potential V is sufficiently large, as one expects 
to occur in the recent epoch. 

In order to simplify the situation, the dynami- 
cal DE is often replaced with models with constant 
w 7^ — 1. This can be considered as a formal gener- 
alization of the equation of state of vacuum energy 
density for which w = — 1. These models result in 
accelerated expansion if w exceeds w —1/3. The 
main advantage of constant w is to yield models 
easier to deal with than the dynamical DE. Al- 
though finding a physical justification for models 
with constant w ^ — 1 is more difficult than for 
the cosmological constant (see, however, Caldwell 
2002), these models are still useful as toy models, 
allowing one to inspect the effects of an accelera- 
tion which is slower than with the vacuum energy. 

In this paper we show how complications with 
the dynamical DE can be overcome if one uses 
suitable approximations, that we provide. Besides 
of allowing an easier treatment of the dynamical 
DE, these expressions also allow us to compare the 
dynamical DE with the models with constant w. 
One of the results of this comparison is that differ- 
ences between constant-u; and dynamical DE are 
significant, being comparable with those between 
ACDM and constant-w. 

The results given in this paper are based on 
a modified version of the CMBFAST code. The 
modifications include effects due to the change in 
the rate of the expansion of the Universe and fluc- 
tuations of the scalar field. Although these fluctu- 



ations rapidly fade, soon after their enter the hori- 
zon, their effect on cosmic microwave background 
anisotropics and polarization is quite significant, 
while they also cause (smaller) modifications of 
the transfer function on large scales. 

In addition, we also estimate the growth of lin- 
ear and non-linear fluctuations of non-relativistic 
matter only. Previously our algorithms were used 
by Mainini, Maccio, & Bonometto (2003, MMB03 
hereafter). 

Making use of these algorithms, in this paper we 
work out: (i) Analytical approximations of the de- 
pendence of the matter density parameter il m on 
the redshift z; (ii) Modifications to run N— body 
simulations of clustering of dynamical DE mod- 
els; (iii) Analytical approximations for the virial 
density contrast A c at any redshift z. Expressions 
derived from the linear theory can also be used 
to compare the observables deduced for dynami- 
cal DE and for constant w. We argue that these 
approximations make an analysis of the dynamical 
DE as simple as for models with constant w. 

2. The virial density contrast 

We start with finding the evolution of the den- 
sity contrast in the top-hat approximation for 
models with DE. Considering a spherical fluctua- 
tion greatly simplifies the analytical and numerical 
treatment of the non-linear problem. Much work 
has been done in this line, starting with Gunn 
& Gott (1972), Gott & Rees (1975) and Peebles 
(1980), who studied the spherical collapse in stan- 
dard CDM (SCDM) models. Lahav et al. (1991), 
Eke et al. (1996), Brian & Norman (1998) and 
others generalized the results to the case of the 
ACDM. If the initial density contrast of a spheri- 
cal perturbation is Aj = 1+Si and its initial radius 
is Ri, then the radius of the perturbation R — rRi 
at later times can be found using the equation: 



= -H? 



2r 3 



r.i 1 
V a 
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(1 + iw)p DE 



2pc 



(4) 



where all quantities with subscript i refer to the 
initial time. In particular fi TOl , and Q, r ^ are the 
density parameters for non- and relativistic matter 
at that time. After slowing down relative to the 
scale factor a(t), the perturbation eventually stops 
at moment t ta , when its radius is R ta - The radius 
R formally goes to zero at <~ 2t ta corresponding 



2 



to redshift z co i. The value of z co i depends on the 
amplitude of the initial fluctuation Si. Instead of 
Si it is however convenient to use the amplitude 
S c as estimated by the linear theory at z co i. For 
SCDM the value of this density contrast is 



1.68 



(5) 



and does not depend on z co i (see, e.g., Coles & 
Lucchin 1995). For other models S c does depend 
on z co i ■ 

In the contraction stages fluctuations heat up 
and, unless kinetic energy can be succesfully radi- 
ated away, contraction will stop when virial equi- 
librium is attained and its size is R v . Requiring 
energy conservation and virial equilibrium we ob- 
tain the following algebric cubic equation 



^3 _ l + y(ata) 

2y{a co i) 

where x = R v /R ta and 



x + 



1 



4y(a co; ) 



= 0, 



(6) 



y(a) = 



1 - Q, m {a) ( Rta 
Aitt m (a) 



Note that the actual radius of the final virialized 
halo is often larger than R v , owing to deviations 
from spherical growth in the real world (Maccio, 
Murantc & Bonometto, 2003). Still, R v is a good 
starting point for statistical analysis. Multiplying 
eq. (6) by 2y and taking then y = (i.e. TO = 1: 
SCDM), we recover that x = 1/2. In general, the 
root x lays slightly below this value. 

Figure 1 shows the linear and non-linear growth 
of a density contrast, for SCDM, ACDM and a 
RP models normalized to have z co i = 0. Similar 
plots can be made for any redshift of collapse. The 
Figure can be used to find the initial amplitude Aj 
at any given redshift Zi and the value of S c for a 
perturbation collapsing at present. Using the final 
value of A we obtain the virial density contrast: 



A, = n m A. 



(8) 



In the linear and non-linear cases deviations 
from the SCDM behavior often compensate and 
the final values of 5 C arc just slightly model depen- 
dent (see Figure 2 and MMB03 for more details). 
The spread among the virial density contrasts A c , 
is large as indicated by Figure 3, which shows A c 
as the function of f2 m for different models. The 




Fig. 1. — Normalized linear (bottom curves) and non- 
linear (top curves) amplitude of density fluctuations 
for SCDM (dotted), ACDM(dashed), and RP (full) 
models. The amplitude of fluctuation was normalized 
to have collapse if the perturbation at z co i — 0. Similar 
plots can be given for collapse at any other redshift. 
The density contrast A = A c /Q m . 



evolution of A c with redshift is also very model 
dependent, as shown by Figure 4). We provide an 
approximation, which is valid at any redshift z, 
provided that we know the matter density param- 
eter Vl m at that redshift: 



A c = 178^ ( ° m ' A) . 



(9) 



Here fj,(Sl m , A) = a + bil^ with c = 1 (2) for RP 
(SUGRA). Parameters a and b are given by 



a = ai\ + a 2 , b = bi\ + b 2 , (10) 



where 



A = log(A/GeV) (11) 

and the coefficients are given in Table 1. 

Figure 5 shows the dependence on A of the dif- 
ferences \A™ m /A% n -l\, at z = 0, for models with 
h = 0.7 and different values of A, as a function of 
TO . (Here A™ um is obtained from the full numer- 
ical treatment, while A"™ is the expression (9)) 
Discrepancies stay below 0.5 % for any fl m ^> 0.15. 
However, for large A, the approximation is even 
better: ^ 0.2%, for any Q m . 
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Fig. 2. — The dependence of 8 C on the matter density 
parameter Q m at z = for 4 RP (A/GcV = 10 2 , 10 4 , 
10 6 and 10 8 ) and 2 SUGRA models (A/GeV = 10 2 and 
10 8 ). A values increase from top to bottom curves. 



Fig. 4. — Redshift dependence of A c for different cos- 
mologies. RP and SUGRA models, at z = 0, have a 
pressure/density ratio w = w of the constant w mod- 
els shown. Full curves are for DM and ACDM. 




Q m 

Fig. 3. — fl m dependence of A c for different cosmolo- 
gies. RP and SUGRA models, at z = 0, have a pres- 
sure/density ratio w = w of the constant w models 
shown. The full curve is for ACDM. 

3. The mass function and the linear growth 
factors for dynamical DE 

We use both the Press-Schechter (PS; 1974) 
and the Sheth-Tormcn (ST; 1999, 2002) approx- 



imations for the mass function of dark matter 
halos. The value of S c defines the bias factor 
v = 5 c /<Jm for the mass M. Here gm is the rms 
density fluctuation on this scale. The bias factor 
enters then the expression 

!{v)vd\ogv = — n h (M)Md log M, (12) 

Prn 

with either 

f(u) v = y/2/w ^exp(-^7 2 ) , (PS) (13) 

or 

f{v)v = 

A(l + !/- 2 <?)y2^rVexp(-z/ 2 /2), (ST) (14) 

with a small complication in the ST case; here 
v' = y/av with a = 0.707, while the constants 



Table 1: Interpolation coefficients for A c 



Model ai 


a 2 


h b 2 


RP -1.45 xlO" 2 
SUGRA -2.25 xl0~ 3 


0.186 
0.3545 


-0.011 0.22 
-0.01875 -0.1225 
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Fig. 5. — Fractional discrepancy between numerical 
and analytical results on A c . 
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Fig. 6. — Linear growth factor for various models. w a 
is the value of w at z — 0. 



q = 0.3 and A = 0.3222. Using cq. (12) we obtain 
the differential mass function (MF) rih(M) in the 
PS and ST approximations, once the distribution 
on bias is given. Here, as usual, we assume a Gaus- 
sian f{v). Eqs. (12-14) can then be integrated to 
obtain the halo mass function nh(> M,z) at any 
redshift z. 

Such computation must use appropriate values 
for 5 C and ctm; the latter are computed by inte- 
grating the power spectrum P(k). Its shape de- 
pends on specific choice of dynamical DE as our 
modified CMBFAST program shows. Yet, the de- 
pendence is very mild for wavelength smaller than 
the galaxy cluster scale. On the contrary, as can 
be seen also from Figure 1, the linear growth fac- 
tor depends on DE nature in quite a significant 
way. Figure 6 presents the z-dependence of the 
growth factor for z up to 40 and for a number of 
different models. 

In particular, Figure 6 shows that at redshift 
z > 2 the difference between ACDM and a model 
with constant w = —0.86 is equal or even smaller 
than the very difference between this constant-w; 
model and the SUGRA model yielding the same 
ratio w at z = 0. However, the latter difference 
becomes comparable with the former one already 
at z ^ 0.5. A constant-w approximation seems 



to perform better for RP models, but this can be 
mostly ascribed to the fact that the ratio w at 
z = is smaller, for these models. Their dis- 
tance from ACDM is therefore greater and the 
difference between them and constant-w models 
appears comparatively smaller. However, also in 
this case, using constant w instead of RP, at z > 4 
is surely misleading. 

The linear growth factor shown in Figure 6 
is very important for setting initial conditions of 
A— body simulations because linear growth fac- 
tors are required for the normalization of the 
power spectrum at the initial redshift z- m of sim- 
ulations. For these reasons we also give an ana- 
lytical approximation which reproduces fairly well 
the behavior of the linear growth factors at z = 40 
for different values of f2 m (z = 0) and A: 



S(z = 40) 



= A + BX + CX 2 . 



(15) 



The values of the coefficients A, B, C are presented 
in Table 2 for RP and SUGRA models respec- 
tively At z — 40, the discrepancies between Q m 
and unity already range around 2-3 %. If a sim- 
ulation must be started at larger z, extrapolating 
the linear growth factor by assuming that S oc a 
at z > 40, implies an error smaller than such per- 
centage. This can be still improved by assuming 
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that 5 oc a 57 ™ , with q ~ 0.4. The dependence of q 
on the model and on the energy scale A fixes the 
second decimal of q and allows a precision bet- 
ter than 0.01 %, which is out of the scopes of this 
analysis. 

Mass functions n(> M,z), obtained according 
to eqs. (12-14), will be compared with simula- 
tions in the accompanied paper (Klypin, Maccio, 
Mainini & Bonometto 2003). Similar mass func- 
tions, obtained from the PS expressions, were used 
in MMB03 to estimate expected observable dif- 
ferences between models with different dynamical 
DE. 

4. Evolution of the matter density param- 
eter 

In RP and SUGRA models, at variance from 
models with w = const, no analytical expression 
of f2 m (a) is readily available. An accurate approx- 
imate expression of Q m , for various redshifts and 
for different models is useful for various purposes. 
In particular, it can be used, in conjunction with 
eqs. (9-10), to find the value of A c at z ^= 0. 

We found the following fitting formula: 

Q m (a) = 1 - (1 - O m , )/(1 + z) tt(z ' A \ (16) 

where fl m ,o is the matter density parameter at z = 
0, while a(z, A) = a + bz c + d/(l + z) with d = 
for RP models. Parameters a, b, c and d have 
the same structure as eq. 10. The coefficients are 
given in Tables 3. 

Figure 7 shows the errors of approximation 
|fi™ m /fi™ - 1| as a function of the rcdshift z 
for two RP and for two SUGRA models with 
Q m = 0.3 and h = 0.7. 



Table 2: Coefficients for the linear growth factor 



Parameter 


n m = 0.2 


On - 0.3 


On = 0.4 


SUGRA 








A 


25.6 


28.5 


30.7 


B 


-0.237 


-0.26 


-0.274 


C 











RP 








A 


21.3 


25.1 


28.2 


B 


-0.755 


-0.783 


-0.698 


C 


-0.0125 


-0.0155 


-0.0155 




Fig. 7. — Fractional discrepancies between the ap- 
proximated expression (16) and numerical data. 

All that is needed to find the relation between 
the scale factor a and time in any flat dynami- 
cal DE model, is such expression. In fact, let us 
remind that 



a/a = H ^Jp(a)/p cr ,o (17) 

with 

p(a) = Pm,o/a 3 + Pr/a 4 + <p 2 /2 + V{<j>) . (18) 

At low z, we can omit the contribution of the ra- 
diation density. Therefore, at any time, 



P4> = 



+ V(4>)= Pcr (a)[l-n m (a)] 



Pm{a) 



1 - fi m (a) 
fi m (a) 



(19) 



provided that we are dealing with a model, such 
that the total density is equal to the critical den- 
sity p cr (a). Then, the Friedmann equation reads 



(") : 
a 



so that 



fl m (a) 



3 a 3 n m (a) ' 




a 3 f2 m (a) 



(20) 



(21) 
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Table 3 
Coefficients for Sl m (z) 



Parameter 


ft m = 0.2 




Mm = 0.3 




tt m = 0.4 




RP 














on 


-5.638 x 10" 


-0 


-2.119 x 10" 


_2 


-3.365 x 10" 


_2 


a 2 


-0.813 




-0.259 




0.207 






-2.460 x 10" 


-2 


-1.833 x 10" 


- 2 


-1.384 x 10" 


-z 




1.382 




0.975 




0.628 




Cl 


-5.960 x 10" 


3 


-6.975 x 10" 


-3 


-8.394 x 10" 


3 


C2 


8.460 x 10~ 2 




9.771 x 10- 2 




0.119 




SUGRA 














Oi 


-8.466 x 10" 


3 


-9.161 x 10" 


-3 


-2.035 x 10" 


-2 


a2 


1.383 




1.415 




1.427 




61 


-1.386 x 10" 


-2 


-1.753 x 10" 


-2 


-1.336 x 10" 


-2 


&2 


-8.521 x 10" 


3 


-6.890 x 10" 


-3 


-1.289 x 10" 


2 


Cl 


-3.935 x 10" 


■2 


-4.421 x 10" 


-2 


-4.203 x 10" 


-2 


C2 


0.710 




0.688 




0.682 




di 


2.088 x 10~ 2 




1.875 x 10~ 2 




2.212 x 10~ 2 




d 2 


-0.883 




-0.621 




-0.416 





This formula is valid regardless of the equation 
of state of DE. In models with constant to, the 
density pde a~^ 1+w ^ and, therefore, owing to 
cq. (21), 

n m (a) = [l + a- 3w (n m 1 -l)}- 1 . (22) 

The expressions (16-22), yielding £l m (a), as 
well as eq. (21), yielding p cr , can be used in 
N— body programs, to determine the trajectories 
of particles in an expanding universe. In fact, once 
we know TO and p cr , we can integrate the Pois- 
son equation V 2 $ = — iirGa 2 p cr Q m 5 m , yielding 
the peculiar potential $ due to the density fluctu- 
ations S m , obtained from the particle distribution. 
Then, the equations of motion of each particle 

dp . „ , dx p . , 

/ = -aV x $, T = -rr 23 

da da a z a 

(see, e.g., Peebles 1980; here p = av) can be inte- 
grated, using a given by eq. (21), and we obtain 
the evolution of particle positions, as a function 
of the scale factor a. The N— body code ART 
(Kravtsov et al. 1997), used in the accompany- 
ing paper (Klypin, Macco, Mainini & Bonometto 



2003) to discuss the evolution of models with the 
dynamical DE and DE with constant w, has been 
modified on these bases. 

Figure 8 compares the expansion law a apx (t), 
obtained using eqs. (16) and (17), with the nu- 
merical behavior a num (t). Discrepancies seldom 
exceed 0.4% and mostly are well below 0.1 %. For 
any practical reasons the errors are negligible. 

5. Discussion 

Observational effects of DE have been consid- 
ered by various authors, but often models with 
a constant w are used. Besides of being sim- 
pler, constant to models give a feeling that re- 
sults are generic in a sense that they do not de- 
pend on the nature of underlying dark energy. For 
instance, Wang & Steinhardt (1998), Steinhardt, 
Zlatev & Wang (1999), Zlatev, Wang & Steinhardt 
(1999) and Lokas (2002) derived the A c depen- 
dence on Q m and to, in the constant to approxima- 
tion. Schuecker et al (2003) extended the results 
to large negative to values to include the case of 
phantom energy (Caldwell 2002, Schulz & White 
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-0.2 -0.4 -0.6 -0.8 -1 



Fig. 8. — Fractional discrepancies between the ana- 
lytical and numerical integration of eq. (17) to obtain 
o(t) 

2001) . 

Unfortunately, results depend on what is as- 
sumed for the DE. Figure 3 shows the dependence 
of A c on w for models with Q m = 0.3 and h = 0.7 
for three cases: DE is cosmological constant, con- 
stant w ^ — 1, and for dynamical DE with RP or 
SUGRA potentials. The difference between con- 
stant w and dynamical DE is as large as the differ- 
ence between ACDM and a constant w. In other 
words, if we need to consider models more sophis- 
ticated than ACDM, it is not enough to discuss 
only constant w ^ — 1. Figure 9 illustrates that 
the growth factor for dynamical DE cannot be ap- 
proximated by a model with constant w. It seems 
clear that the Universe "knows" the underlying 
physics, and predictions depend on the shape of 
the potential of the scalar field responsible for the 
DE. 

In principle, finding astrophysical quantities 
which depend on microphysics is far from being 
unwelcome. Accordingly, the detailed dependence 
of astrophysical observables on microphysical pa- 
rameters deserves to be inspected. Let us outline, 
in particular, that Figure 3 applies to observations 
at z = 0, while effects of dynamical DE are also 
expected at higher z. In fact, in Figure 4 we show 
the z dependence of A c , for 3 sets of models char- 



Fig. 9. — The linear growth factor for models with 
w = const (dotted curve) is compared with linear 
wrowth factor for RP and SUGRA models (dashed and 
dot-dashed curves, respectively). RP and SUGRA re- 
sults are plotted as a function of the value of w they 
have either at z = or at z — 40 (long and short 
dashes, respectively). The logarithmic energy scale 
A, for both models, ranges here from 2 to 10. The 
plot illustrates that the growth factor for dynamical 
DE cannot be easily approximated by any model with 
constant w. 

acterized by the same values of w at z — 0. The 
figure shows that the differences between A c 's in- 
crease from z — toward intermediate redshifts, 
to go back to SCDM values at high redshifts, when 
ordinary matter gradually approaches critical den- 
sity. However, at intermediate redshifts, which are 
most relevant for present and future observations 
and for the actual dependence on the nature of 
DE, at these z, arises from actual changes of w, as 
shown in Fig. 9. Apart of any consideration con- 
cerning fundamental physics, therefore, high-z A c 
values obtained within constant w approximation 
risk to create a bias. 

Our aim is to facilitate the usage of dynamical 
DE. We provide the following tools: (i) An ap- 
proximation for f2 m (o); (ii) An interpolating ex- 
pression for A c , valid at any redshift for given 
£l m (a); (iii) An analytical expression for the rate 
of change of the expansion parameter needed for 
running N— body and hydro- simulations; (iv) A 
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Fig. 10. — Redshift dependence of w for 4 RP and 4 
SUGRA models (A = 1, 3, 5 and 7); A decreases from 
top to bottom curves. 

plot of the linear growth factor, for a number of 
dynamical DE models, and an analytical approx- 
imation for it, to be used to set the initial condi- 
tions of N— body simulations. 

Using these formulae, we modified the ART 
code, that will be used in the accompanying pa- 
per (Klypin, Maccio, Mainini & Bonometto 2003). 
In a similar way, other programs dealing with 
A^— body interactions or hydrodynamics can be 
appropriately modified. 
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